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ABSTRACT 


This  thesis  presents  a  computer  program  designed  to 
analyze  the  data  from  a  TV  camera  for  investigating  laser 
beam  propagation  through  the  atmosphere.  It  uses  aspects  of 
Fourier  optical  theory  to  analyze  the  TV  image  to  measure 
the  effects  of  atmospheric  disturbance  and  platform  stab¬ 
ility  on  the  target  spot.  The  understanding  and  analysis  of 
these  effects  are  increasingly  important  as  lasers  and  other 
optical  devices  experience  increased  use  in  military 
applications . 
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I.  INTRODUCTION 


A.  PROBLEM 

A  problem  has  occasionally  been  experienced  in  concen¬ 
trating  laser  energy  from  a  laser  designator  onto  a  target 
and  sufficiently  illuminating  that  target  so  a  weapon 
sensitive  to  that  illumination  will  home  in  on  the  target. 
Two  circumstances  have  possibly  caused  this  to  happen;  the 
energy  is  either  too  diffused  at  the  target  to  illuminate  it 
properly  or  the  beam  partially  or  completely  misses  the 
target.  The  aim  of  this  study  is  to  develop  a  computer 
program  to  aid  in  determining  which  part  of  these  effects  is 
atmospheric  and  which  is  the  instability  of  the  laser  and 
stabilization  system  in  the  illuminating  aircraft. 

B.  METHOD  OF  APPROACH 

The  analysis  of  the  problem  requires  two  general  types 

of  systems.  One  system  evaluates  the  performance  of  the 

designator  using  some  techniques  from  Scott  [Ref.  1].  A 

silicon  vidicon  views  the  successive  pulsed  designator  spots 

on  the  target  from  a  location  near  the  target.  A  line 

spread  function  of  the  spot  on  the  target  and  the  standard 

deviation  of  the  wander  of  the  laser  beam  is  produced.  The 

2 

other  system  provides  measurement  of  Cn  for  atmospheric 
turbulence  along  an  optical  path  adjoining  the  designator 


5 


V 


optical  path.  The  laser  on  the  target  is  viewed  by  a 
silicon  vidicon  located  near  the  designator’s  optics.  The 

p 

value  of  Cn  obtained  is  inherently  properly  path  weighted 

to  express  the  atmospheric  effects  on  the  designator  spot. 

Techniques  from  Fourier  optics  theory  use  the  measured  value 
2 

of  Cn  to  predict  the  intensity  pattern  of  the  laser  beam  on 
a  distant  target.  The  results  from  the  two  systems  can  then 
be  analytically  compared  to  determine  the  causes  of  spot 
wander  and  broadening. 

The  sponsor-supplied  hardware  includes  an  HP-9825 
calculator  with  23K  bytes  of  internal  memory  for  equipment 
control  and  data  processing;  an  HP-9885  disc  memory  used  for 
additional  storage  space;  and  a  Quantex  DS-30  Digital  Video 
Analyzer  used  for  digitizing  the  analog  video  input  from  an 
Eigen  video  disc;  a  Panasonic  NV-1240  video  tape  recorder 
for  initial  recording  of  video  from  the  target  on  a  silicon 
vidicon  camera.  After  analysis,  output  is  in  the  form  of 
plots  produced  on  an  HP-9862  plotter. 

The  sequence  of  analysis  is  currently  as  follows. 
Modulation  transfer  functions  (MTF)  of  the  laser  output  and 
optics  are  measured  or  calculated  and  stored  on  disc.  A 
video  recording  of  the  target  being  illuminated  by  the 
designator  is  made  for  a  number  of  video  frames,  the  target 
alone  is  also  recorded,  for  later  subtraction  from  the  total 
picture  to  produce  an  image  of  the  laser  spot  alone.  The 
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video  recording  is  then  played  through  the  DS-30,  controlled 
by  the  HP-9825,  to  digitize  the  information  for  use  by  the 
HP-9825.  The  HP-9825  takes  the  digitized  video  of  the 
target,  produces  a  line  spread  function  (LSF)  of  the  image, 
subtracts  the  background,  and  averages  a  number  of  frames  to 
produce  a  short  term  measurement  of  the  laser  energy  on  the 
target.  The  amount  that  each  laser  spot  wanders  from  the 
others  is  also  stored.  In  the  second  phase,  the  HP-9825 
uses  the  MTF's  of  the  laser,  optics,  and  atmospheric 
measurements  to  produce  a  total  system  predicted  MTF,  LSF, 
and  wander.  These  results  of  the  target  spot  measurement 
and  the  baseline  studies  can  be  compared  for  possible 
correlation  . 

The  laser  pulse  is  timed  so  that  it  occurs  on  the  fly¬ 
back  of  the  TV  signal.  This  will  cause  the  laser  spot  to 
appear  on  one  field  of  a  frame.  For  accurate  results  the 
DS-30  needs  to  digitize  a  single  field  for  each  pulse  of  the 
laser.  The  videotape  unit  currently  used  will  not  allow 
this  capability.  For  this  reason  the  tape  recording  video 
is  transferred  to  an  Eigen  video  disc,  which  has  the 
capability  of  displaying  each  frame  or  field  (1/2  video 
frame  interlace)  individually  and  of  selecting  any  one  frame 
or  field  for  viewing. 
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II .  BACKGROUND  AND  PROGRAM  DESCRIPTION 

A.  THEORETICAL  DESCRIPTION 

In  order  to  combine  the  individual  optical  components  of 
the  laser  to  the  target  system,  some  elements  of  Fourier 
theory  must  first  be  remembered.  The  Fourier  transform 
represents  the  one  dimensional  position  variable  g(y) 
expressed  in  the  spatial  frequency  domain  U(v).  The  inverse 
transform  repeats  the  operation  in  the  opposite  direction. 
Their  forms  are  commonly  represented  as  follows: 


f  g(y)exp(-2 

t  ivy)dys 

g(y)=  1 

f  U(v)exp(2  ir 

ivy) dv  = 

.TCgCy)] 

1  [U  ( v )  ] 


The  program  calculates  the  value  of  the  integral  at  a  preset 
number  of  points  and  yields  a  discrete  Fourier  transform.  A 
theorem  from  Fourier  theory  that  allows  us  to  compute  the 
line  spread  function  as  presented  in  Fried  [Ref.  2]  for  one 
of  the  elements  in  the  system  given  the  others  is  the  con¬ 
volution  theorem  stated  as  follows: 

If  «*■  C g(y )  ]  =  G(v) 
and  dF  [h(y) ]  =  H( v  ) 
then  &  [g*h]  =  GH 

(*  -  convoluted  with) 
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The  first  element  in  the  system  we  are  analyzing  is,  of 
course,  the  laser  itself.  The  radial  intensity  output  of 
the  laser  is  generally  Gaussian  of  the  form  A  =  A  expC-r^/^cj  ) 
as  shown  in  Figure  1.  This  can  be  calculated  for  each  laser 
for  approximate  results  but  was  usually  measured  directly  in 
this  study. 

The  above  radial  distribution  of  intensity  must  then  be 
integrated  over  x  to  produce  a  one  dimensional  line  spread 
function  for  the  source.  A  Fourier  transform  is  then  used 
on  this  LSF  to  yield  the  MTF  of  the  laser.  The  MTF  is  a 
measure  of  the  spatial  frequency  response  of  a  system 
compared  to  the  input.  A  "perfect"  system  response  would  be 
1.0  out  to  its  limit  (large  spatial  frequencies)  and  then 
dropping  to  zero. 

The  optics  associated  with  the  laser  is  basically 
diffraction  limited.  The  "Airy  function"  is  used  to  cal¬ 
culate  the  diffraction  point  spread  function  for  the  laser 
optics.  This  point  spread  function  is  then  converted  to  a 
line  spread  function,  Fourier  transformed  to  get  the  MTF  and 
then  combined  by  the  convolution  theorem  into  a  laser  system 
Fourier  transform. 

The  next  step  is  to  calculate  the  effects  of  atmospheric 
turbulence  on  the  laser  beam.  The  value  of  C^,  "the  index 
of  refraction  turbulence  structure  constant,"  as  expressed 
by  the  relationship  from  Tatarski  [Ref.  33  and  Ochs  et.  al. 
[Ref.  43: 
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FIGURE  1.  IDEAL  LASER  INTENSITY  OUTPUT  LINE  SPREAD  FUNCTION 


<(Trl"Tr2)2  >=r2/3ci2  and  Cn=79x 1 0‘6p/T2CT 

where  I\ , -T  _  is  absolute  temperature  at  two  points 
separated  by  a  distance  r  and  p  is  pressure  in  millibars;  is 
used  to  predict  the  MTF  of  the  atmosphere  over  the  path 
between  the  laser  and  the  target.  Measurements  made  along 
the  path  are  inputs  to  the  Cn  equation  which  predicts  the 
turbulence  effects  on  the  MTF  and  wander  of  the  beam.  In 
the  field,  however,  plans  are  to  measure  Cn  by  optical 
methods  decribed  in  Crittenden  [Ref.  5]. 

The  theory  then,  according  to  Fried  [Ref.  2],  shows  that 
the  total  effect  of  the  atmosphere  and  the  laser  can  be 
found  by  the  product  of  the  two  transform  functions: 

G  total  (v)=G  optics(v)  x  G  atm(v) 

The  inverse  transform  then  will  yield  the  line  spread 
function  of  the  total  system,  the  LSF  as  predicted  on  the 
target . 

Crittenden  et .  al ,  [Ref.  5]  describe  the  process  of 
converting  the  one  dimensional  LSF  to  a  circular  symmetric 
PSF  using  the  Abel  transform,  from  Griem  [Ref.  6].  The 
power  inside  a  circle  of  radius  R  can  then  be  obtained  by 
integrating  the  PSF  out  to  R.  For  the  purposes  of  compari¬ 
son  in  this  study,  however,  the  resulting  system  LSF  is  used 
to  compare  the  prediction  with  measured  values. 
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B.  HARDWARE  USE  AND  LIMITATIONS 

The  HP-9825  calculator  used  for  this  program  had  an 
internal  capacity  of  23K  bytes.  This  was  quite  adequate  for 
program  storage  but  because  of  the  great  volume  of  digitized 
video  data  that  was  involved,  the  inclusion  of  an  HP-9885 
disc  drive  system  was  necessary.  The  Quantex  DS-30  can 
store  and  rapidly  transfer  a  maximum  of  512  lines  of  video 
with  512  pixels  (picture  elements)  per  line,  easily 
exceeding  the  HP-9825’s  memory  capability.  The  large 
overflow  of  data  was  reduced  by  storing  each  image  into  a 
quarter  of  the  memory,  thus  limiting  the  amount  of  data  that 
needed  to  be  transferred.  This  also  made  disc  storage  more 
reasonable  and  sped  up  program  running  time. 

The  DS-30  was  also  capable  of  taking  the  difference  on  a 
pixel  by  pixel  basis  between  a  reference  image  and  the  input 
and  storing  that  in  memory.  Differencing  was  necessary  in 
order  to  remove  the  background  from  the  recorded  image  for 
processing  the  laser  spot  alone.  Normally,  it  is  necessary 
to  divide  the  image  by  the  background  to  offset  the  effect 
of  the  non-uniform  reflectivity  of  the  target.  This  was  not 
possible  with  the  DS-30  so  these  steps  need  to  be  performed 
by  the  computer.  The  division  was  not  possible  to  do  at 
this  time  because  of  the  memory  limitations  of  the  HP-9825. 
This  must  be  done  pixel  by  pixel  since  division  by  the  LSF 
is  not  equivalent.  The  memory  limitation  was  also  the 


reason  why  the  HP-9825  used  the  LSF  instead  of  pixel  by 
pixel  subtraction  for  this  process. 

All  the  peripheral  equipment  to  the  HP-9825  was 
controlled  through  an  IEEE  standard  8  bit  interface  bus. 

This  allowed  control  of  all  aspects  of  data  acquisition  and 
processing  to  be  modified  by  software.  The  video  recorders 
were  the  exception  to  this  and  provided  the  only  real  manual 
manipulations  required  after  the  program  had  begun.  Each 
frame  of  video  needed  by  the  program  had  to  be  selected  when 
requested  by  the  program  because  of  the  lack  of  interface 
control . 

The  linearity  of  the  recording  system  was  checked  and 
found  to  be  almost  distortion  free.  A  signal  from  an  image 
was  passed  directly  to  an  oscilloscope  from  the  vidicon  and 
compared  side  by  side  to  a  signal  from  the  same  image  which 
had  been  recorded  in  the  tape  recorder,  transferred  to  the 
video  disc  and  then  passed  through  the  DS-30  circuitry  to  a 
digital  analog  converter  for  viewing.  Crittenden  et .  al. 
[Ref.  7]  present  photographs  of  the  two  images  for 
comparison  and  further  discussion. 

Variable  usage  became  a  problem  as  the  program  grew, 
which  was  due  to  the  limitations  of  the  HP-9825.  The 
HP-9825  has  available  twenty-six  variables  (A-Z),  twenty-six 
variables  for  arrays  (A-Z),  a  subscripted  variable  (r)  with 
as  many  elements  as  the  memory  has  space  for  plus  a 
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subprogram  variable  (p)  for  use  in  passing  parameters 
through  "called'*  subprograms.  One  of  the  usage  problems  was 
using  identifiable  variables  in  passing  from  main  programs 
to  subprograms.  A  good  number  of  the  subprograms  originally 
were  designed  for  other  programs  and  the  variables  often 
conflicted.  The  r  variable  is  nice  to  use  here  because  of 
its  practically  unlimited  numbers.  A  more  serious  problem, 
though,  was  the  inability  of  HPL,  Hewlett  Packard's  version 
of  the  Basic  language,  to  allow  an  entire  array  to  be  used 
as  a  parameter  to  pass  to  a  subprogram.  Thi3  at  first  made 
necessary  extra  internal  storage  for  additional  arrays  to  be 
used  only  within  a  certain  subprogram,  with  the  variable 
name  being  changed  through  a  lengthy  process  upon  entering 
or  leaving  that  subprogram.  A  much  faster  and  easier  way 
was  using  separate  storage  files  on  the  disc  for  each  major 
array  manipulation.  When  entry  into  a  subprogram  was  needed 
the  array  could  be  stored  in  the  file  assigned  to  that 
subprogram,  which  could  then  be  read  out  and  labelled  as  any 
array  needed  for  each  subprogram. 

The  programming  to  use  the  DS-30  also  was  primarily  a 
language  problem.  The  DS-30  was  very  sensitive  to  the 
format  in  which  it  received  data  or  commands.  All  leading 
zeroes  and  carriage  return/line  feed  commands  had  to  be 
suppressed.  The  address  of  command  locations  within  the 
DS-30  also  had  to  be  encoded  in  hexadecimal  nibbles  (one- 


half  byte)  and  sent  in  reverse  order  to  load  the  registers 
properly.  Commands  to  the  DS-30  were  sent  as  ASCII 
character  strings.  Data  was  output  from  the  DS-30  to  the 
HP-9825  via  a  fast  read  buffer  which  allowed  a  very  fast 
data  transfer  rate.  In  use,  a  block  of  256  bytes  at  a  time 
were  transferred  representing  one  horizontal  TV  line  of  the 
quartered  memory  space.  This  was  done  in  a  loop  which  went 
back  and  took  256  lines,  completely  transferring  the  image. 

It  was  decided  that  output  would  be  on  the  HP-9825 
plotter  so  that  direct  visual  comparisons  could  be  made 
between  measured  and  predicted  intensity  LSF’s.  This  was 
much  more  dramatic  and  easier  to  see  than  tabular  output 
which  then  would  have  to  be  compared  on  a  point  by  point 
basis  for  a  correlation. 
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III.  CONCLUSION 


This  program  can  be  useful  in  solving  the  problem  that 
spurred  this  study.  It  will  be  able  to  show  how  much  wander 
can  be  expected  out  of  a  laser  designator  beam  as  it  travels 
through  the  atmosphere  from  the  laser  to  the  target.  Based 
on  the  video  recording  of  the  actual  image  on  a  target  it 
will  show  comparatively  the  amount  of  wander  actually 
present  in  the  laser.  Plans  are  to  unfold  the  results 
analytically,  but  at  present,  a  visual  comparison  of  the 
measured  and  predicted  patterns  is  made  to  determine  if 
there  is  any  excess  wander  or  "jitter"  attributable  to  the 
laser  platform  itself.  This  determination  will  be  helpful 
in  setting  specifications  for  the  stability  needed  in  the 
system  to  remove  as  much  excessive  spot  movement  as 
possible.  It  will  also  allow  testing  to  determine  if  those 
specifications  are  being  met. 

This  program  can  also  be  used  to  verify  atmospheric 
problems  when  both  the  laser  and  target  spot  are  fixed  and 
known.  Essentially,  any  of  the  laser  to  target  elements  can 
be  determined  by  knowing  the  others  and  calculating  the 
unknown  using  adaptations  of  the  program.  The  beginning 
part  of  this  program  can  be  used  to  digitize  and  store  any 
video  image  and  produce  an  LSF  of  the  image.  This  may  be 
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useful  in  recording  and  testing  the  output  of  several 
different  lasers  and  optical  aperatures  for  comparative 
analysis . 

This  computer  program  was  designed  in  the  initial  stages 
of  a  continuing  project  studying  the  described  problem.  If 
desired,  changes  can  be  made  to  the  program  to  adapt  to 
changing  field  conditions  and  analysis  requirements.  In 
particular,  the  program  now  requires  a  greater  understanding 
of  its  internal  operations  and  operator  interaction  than  is 
necessarily  needed.  More  complete  cueing  prompts  and  data 
entry  parameters  may  be  worthwhile  changes.  Adaptability 
to  other  hardware  may  provide  more  program  control,  especi¬ 
ally  if  control  of  the  video  disc  or  tape  can  be  handled 
through  an  interface  bus  instead  of  manually.  Also, 
improvements  in  output  design  can  be  foreseen  to  expand  its 
usefulness  and  adaptability. 

Finally,  the  speed  of  program  operation  could  be 
increased  through  use  of  a  faster  computer  with  a  greater 
memory  capacity.  The  HP-9825  is  comparatively  slow  and  has 
insufficient  memory.  The  use  of  a  faster  system  such  as  the 
HP-1000  is  recommended  for  actual  measurements. 
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APPENDIX  A 


USING  THE  COMPUTER  PROGRAM 

The  analysis  program  requires  considerable  preliminary 
set  up  and  operator  interaction.  Before  the  program  is 
started  there  should  be  several  items  already  recorded  on 
videotape.  The  first  item  is  a  TV  camera  shot  of  the  laser 
output.  This  is  normally  taken  with  a  long  focal  length 
lens  and  filters  to  reduce  the  intensity  without  losing  the 
beam  pattern.  A  part  of  this  recording  should  include  a 
known  diffraction  grating  through  which  the  beam  passes  in 
order  to  measure  the  scale  of  the  data  for  each  system  of 
lenses  used.  The  diffraction  pattern  should  extend  in  the 
vertical  as  seen  on  a  TV  screen  for  proper  output  on  the 
plotter.  Also  needed  to  be  recorded  is,  of  course,  a 
sequence  of  laser  spot  images  on  the  target.  Part  of  this 
recording  should  have  the  target  only  as  a  background 
reference  for  later  calculations. 

The  next  step  that  needs  to  be  taken  is  to  compute  and 

2 

record  the  square  atmospheric  turbulence  index  (Cn  ),  then 
record  the  following  system  data  for  later  program  entry; 
wavelength  of  the  laser,  diameter  of  the  optics  objective 
lens,  ratio  of  the  obscuration  to  the  diameter  of  the 
objective  lens,  distance  to  the  target,  and  the  extinction 
coefficient . 
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The  first  cue  the  program  asks  is  for  the  input/output 
parameters.  Table  1  lists  the  parameters  and  the  different 
function  each  one  is  used  for. 

If  the  laser  source  pattern  is  not  to  be  computed  the 
next  cue  requests  that  the  laser  pattern  be  displayed  from 
the  video  disc  for  transfer  to  HP-9825  memory.  A  cue  is 
also  presented  asking  for  the  horizontal  width  of  the  images 
to  be  taken  from  the  DS-30.  The  first  value  should  be  the 
left  most  pixel  desired,  the  second  being  the  furthest  right 
pixel  to  be  used.  The  first  value  must  be  less  than  the 
second,  with  the  first  no  less  than  1  and  the  second  no  more 
than  256.  This  cue  is  also  presented  below  before  either 
the  scale  image  or  the  target  image  is  taken. 

If  the  scale  of  the  source  data  is  to  be  determined  the 
next  cue  requests  that  a  frame  be  displayed  on  the  video 
disc  of  the  laser  output  with  the  diffraction  grating  in 
place.  When  this  has  been  done  continue  the  program  and  it 
will  plot  out  an  LSF  of  that  image  whereby  the  scale  of  the 
data  can  be  determined  in  raicroradians  per  point. 

The  next  cue  alerts  the  operator  to  ensure  that  a  frame 
of  background  video  is  displayed  from  the  video  disc  for  use 
in  the  difference  process. 

The  program  will  then  ask  for  the  number  of  total  frames 
that  the  user  wishes  to  average  for  the  measured  target  spot 
LSF.  The  value  should  be  in  the  range  from  1  to  600.  The 


latter  value  being  the  greatest  number  of  frames  the  video 
disc  can  hold  at  one  time.  The  program  will  then  proceed  to 
record  and  process  images  from  the  video  disc  and  DS-30* 

The  program  will  provide  a  cue  when  it  wants  the  next  image 
to  be  sequenced  on  the  video  disc. 

The  program  will  alert  the  operator  when  it  is  about  to 
produce  a  plot.  If  the  HP-9862  plotter  is  not  ready  simply 
stop  the  program  until  a  piece  of  plotting  paper  is  in 
place . 

The  program  will  request  the  system  data  after  the 
measured  image  spot  data  has  been  plotted.  The  data  format 
is  not  critical  but  should  be  in  the  correct  units  as 
follows : 


Wavelength 

Diameter  of  objective  lens 
Obscuration/objective 
Scale  of  data 
Range  to  target 


Extinction  coefficient 


Meters 

Meters 

None 

Microrads  per  point 
Meters 
Meters-273 
Inverse  meters 


If  the  source  pattern  is  to  be  computed  the  program 
requests  peak  amplitude  and  standard  deviation  of  the 
intensity  distribution . 
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The  program  will  then  continue  computations  without 

further  inputs  until  it  has  run  through  with  the  initial 

parameters.  The  I/O  parameters  will  determine  at  which 

steps  output  is  plotted.  After  one  run  through,  the  program 

2 

will  request  a  new  Cn  ,  range,  and  coefficient.  If  the 
input  range  is  positive  the  program  will  again  run  through 
from  the  atmospheric  calculations  with  the  new  range.  If 
the  range  is  negative,  the  branch  is  to  the  beginning  of  the 
prediction  phase  where  all  new  system  inputs  are  requested. 
If  the  range  input  is  zero,  the  program  ends. 


APPENDIX  B 


PROGRAM  VERIFICATION  AND  FLOWCHARTS 
I.  MEASUREMENT  PHASE 

Section  I  (steps  0  -  50)  consisted  of  program  setup 
information  and  the  DS-30  function  command  programming. 

This  section  was  easily  verified  by  observing  that  all 
initializing  functions  were  accomplished  and  the  response 
of  the  DS-30  could  be  visually  noted. 

Section  II  (steps  51  -  86)  consisted  of  transferring 
the  contents  of  the  four  quarters  of  the  DS-30' s  memory  to 
the  HP-9825  and  producing  an  LSF  of  each  image.  For 
testing,  an  image  consisting  of  all  black  (8  bit  binary 
ones)  or  all  white  (8  bit  binary  zeroes)  was  placed  in  each 
of  the  four  quarters  of  memory,  each  quadrant  was  then  read 
out  to  the  HP-9825.  The  values  received  by  the  HP-9825 
corresponded  to  the  contents  of  the  DS-30  memory  quadrant 
and  did  show  proper  changes  when  the  subject  quadrant 
memory  was  inverted  from  black  to  white.  The  LSF  was 
produced  by  averaging  the  pixels  in  each  horizontal  line  of 
video,  each  line  producing  a  point  value  as  the  image  was 
scanned  vertically  from  top  to  bottom.  Examples  of  LSF's 
produced  for  a  "black"  and  a  "white"  image  are  in  Figures  2 
and  3* 
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Section  III  (steps  87  -  118)  subtracted  a  reference 
background  from  an  image  to  yield  a  difference  image,  then 
shifted  the  resulting  LSF  so  that  the  center  of  area  under 
each  image  was  centered  on  zero.  In  this  way  the  wander 
from  each  succeeding  image  could  be  determined,  removed, 
and  the  total  frames  used  could  then  be  averaged  to  show 
the  short  term  atmospheric  action  on  the  laser  spot  image. 
Proper  operation  of  this  section  can  be  seen  in  Figure  4 
where  two  images  of  slightly  different  position  and 
intensity  were  averaged  to  yield  an  LSF  halfway  between  the 
two  image  functions  as  expected.  Figure  5  shows  a  differ¬ 
enced  LSF  where  a  small  laser  spot  was  added  to  a  reference 
background.  The  result  is  almost  entirely  the  laser  spot 
with  a  small  amount  of  "video  noise"  caused  by  the  camera 
and  associated  electronics. 
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Image 
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FIGURE  4.  AVERAGING  EXAMPLE 


FIGURE  5.  LASER  SPOT  LSF 


rr 
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II.  PREDICTION  PHASE 

Section  IV  (steps  119  -  139)  enters  the  parameters  of 
the  atmosphere  and  optics  needed  to  model  and  predict  the 
laser  spot  pattern  on  the  target.  This  section  also  yields 
the  calculated  Gaussian  laser  intensity  output  if  that  data 
is  not  measured  in  the  next  section.  The  output  of  this 
section  is  shown  in  Figure  7. 

Section  V  (steps  140  -  172)  consists  mainly  of 
branching,  data  reading,  and  subroutine  action.  These  sub¬ 
routines,  branching  and  data  input  will  be  verified  with 
proper  operation  of  the  program  as  a  whole. 

Section  VI  (steps  173  -  183)  calculates  the  diffraction 
limited  optics  effect  on  the  laser  output  and  produces  a 
radial  point  spread  function  of  the  intensity.  The  output 
for  a  point  source  should  be  an  airy  function  as  results 
show  in  Figure  3. 

Section  VII  (steps  184  -  204)  again  consists  mainly  of 
subprograms  and  also  calculates  the  product  of  the  calcu¬ 
lated  optics  and  laser  output  transforms.  The  output  is  in 
Figure  9. 

Section  VIII  (steps  205  -  220)  computes  the  MTF  of  the 
atmosphere  based  on  the  inputs  from  Section  IV.  The  output 
is  in  Figure  10. 

Section  IX  (steps  221  to  the  end  of  the  main  program) 
is  mainly  subprogram  action  in  computing  the  final  LSF 
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predicted  on  the  target  using  the  source,  optics  and 


atmospheric  MTF's.  It  also  provides  a  link  to  re-enter  the 
program  for  changed  parameters.  The  output  of  this  section 
is  listed  under  the  subprograms  used. 
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FIGURE  7.  EXAMPLE  OF  CALCULATED  LASER  OUTPUT 


FIGURE  8.  OPTICS  OUTPUT 


FIGURE  113.  PREDICTION  PHASE  FL' 
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III.  SUBPROGRAMS 

Section  X  (subprogram  ABEL)  produces  a  point  spread 
function  (PSF)  from  the  LSF  computed  for  the  target  laser 
spot.  The  input  and  resulting  output  is  shown  in  Figures 
12  and  13. 

Section  XI  (subprogram  LSF)  enables  the  program  to 
compute  a  line  spread  function  for  later  conversion  to  a 
MTF  when  the  input  is  calculated  by  Sections  IV  and  VI  as 
point  spread  functions  or  when  the  data  is  input  in 
Section  V  in  the  point  spread  function  format.  Sample 
input  and  output  is  in  Figures  14  and  15. 

Section  XII  (subprogram  B)  calculates  the  fraction  of 
the  power  inside  a  circle  of  a  chosen  radius  about  the 
center  of  the  laser  spot.  This  is  helpful  in  determining 
if  the  spot  is  intense  enough  on  the  target.  The  input  and 
output  are  Figures  16  and  17. 

Section  XIII  (subprogram  FXFORM)  is  a  discrete  Fourier 
transform  program  that  functions  both  as  direct  and  inverse 
depending  on  the  flag  set  in  the  main  programs.  Several 
Fourier  pairs  were  checked  to  ensure  correct  calculations, 
Figures  18  through  21. 

Section  XIV  (subprogram  SOURCE)  digitizes  an  image  of 
the  laser  source  for  inclusion  in  the  prediction 
calculations.  It  is  a  storage  of  DS-30  video  from  the 
Vidicon  camera  input. 
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Section  XV  (subprogram  SCALE)  digitizes  an  image  with  a 
standard  diffraction  grating  for  calculating  the  scale  of 
the  data.  Figure  22  is  an  example  of  this  image. 

Section  XVI  (subprogram  OUTPUT)  yields  in  useable  form 
all  the  data  required  for  analysis.  This  subprogram  was 
used  for  all  the  output  in  the  plot  format.  Functions 
"MIN"  and  "MAX"  are  used  by  the  "OUTPUT"  subprogram. 

The  program  was  run  as  a  whole  numerous  times  with 
different  I/O  parameters  and  data  in  order  to  exercise  the 
branching  combinations  and  program  steps.  Table  1  lists 
the  I/O  parameters  available  and  their  uses.  This 
operation  was  simply  to  ensure  the  program  ran  as  was 
intended . 
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FIGURE  22.  SCALE  IMAGE  EXAMPLE 


TABLE  1 


INPUT/OUTPUT  PARAMETERS 


ARRAY  VARIABLE 

VALUES 

USE 

0(1) 

0 

Same  as  1  except  program  after  one 
time  through  ends,  not  allowing2for 
computations  with  new  range,  C  , 
ext.  coeff.  1 

1 

Source  pattern  input  is  point 
spread  function. 

2 

Source  pattern  input  is  line  spread 
function  . 

3 

Source  pattern  input  is  modulation 
transfer  function  of  laser  only. 

4 

Source  pattern  input  is  modulation 
transfer  function  of  source  and 
optics . 

5 

Source  pattern  is  to  be  calculated 
with  amplitude  and  sigma  of  inten¬ 
sity  distribution  as  factors. 

0(2) 

0 

Optics  diffraction  input  as  PSF. 

1 

Optics  diffraction  input  as  LSF . 

2 

Optics  diffraction  pattern  to  be 
calculated  as  PSF. 

0(3) 

0 

No  plot  of  source  line  spread 
function  is  required. 

>0 

Plots  source  line  spread  function 
(available  only  if  0(1 )=0, 1  or  5). 

0(4) 

0 

No  plot  of  source  MTF  required. 

>0 

Plots  source  MTF  (available  only  if 
0(1)  does  not  equal  3  or  4). 

0(5) 

0 

No  plot  of  product  of  source  and 
optics  MTF  required. 

>0 

Plots  product  of  source  and  optics 
MTF's  (available  only  if  0(1)  does 
not  equal  4). 
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TABLE  1  (CONT'D) 


ARRAY  VARIABLE 

VALUE 

USE 

0(6) 

0 

No  plot  of  optics  diffraction  PSF 
required . 

>0 

Plots  optics  diffraction  PSF 
(available  only  if  0(1)  does  not 
equal  4  and  0(2)=2). 

0(7) 

0 

No  plot  of  optics  LSF  required. 

>0 

Plots  optics  LSF  (available  only  if 
0(1)  does  not  equal  4  and  0(2)  does 
not  equal  1 ) . 

0(8) 

0 

No  plot  of  optics  MTF  required. 

>0 

Plots  optics  MTF  (available  only  if 
0(1)  does  not  equal  4). 

0(9) 

0 

No  plot  of  atmosphere  MTF  required. 

>0 

Plots  atmosphere  MTF. 

0(10) 

0 

No  plot  of  system  MTF. 

>0 

Plots  total  system  MTF. 

0(11) 

0 

No  plot  of  LSF  of  predicted  target 
spot . 

>0 

Plots  LSF  of  target  spot. 

0(12) 

0 

No  plot  of  target  PSF  required. 

>0 

Plots  target  PSF. 

0(13) 

0 

No  plot  of  the  fraction  of  power 
inside  circle  required. 

>0 

Plots  fraction  of  power  inside  a 
circle  of  radius  r. 

0(14) 

0 

any 

other 

value 

No  plot  of  calculated  source  PSF 
required . 

Plots  calculated  source  PSF. 
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TABLE  1  CCONT  *D) 


ARRAY  VARIABLE 

VALUES 

USE 

0(15) 

0 

Plots  measured  LSF  of 

target  spot. 

>0 

Used  as  indicators  of 
branching  only,  entry 
significant . 

output 
value  not 

0(16) 

1 

Inputs  source  with  diffraction 
grating  from  DS-30  for  scale  compu¬ 
tations  from  output.  No  output 
plotted . 

2  Same  as  1  except  plots  the  above. 


any 

other  Does  not  do  above  input 
value 


APPENDIX  C 


PROGRAM  LISTING 


0:  dim  B$[23]fC$l5],D$[7]  ,I$[272]  ,0$[5],P$15J  ,si*lS],X$l4Ql  fR[ 

512]  ,1.(512] 

1:  dim  0116]  ,AR[13]  ,¥$.(4Q]  ,ZS(40] 

2:  files  O rT,£ ,R ,G ,V , W,C, I ,? 

.3:  f mt  1,  z,  c 

4:  time  150 

5:  cfg  3;cfg  2;cfg  1 

6:  rread  7,l;sprt  7,R[*];rread  7,1 

7:  buf  "DATA", 1$, 3 

8:  rread  4,l;sprt  4,R[*] 

9  •  p 

10:  "INPUT  I/O  CONTROL  0(1]": 

11:  for  1*1  to  16;ent  0(I];next  I 
12:  if  0 (1 ] ®2 ; gso  "SOURCE" 

13:  clr  704 

14:  if  0  [16]  =1  or  0[16]=2;gsb  "SCALE" 

15:  beep;dsp  "ENSURE  VIDEODISC  ON  THEN  PRESS  CONTINUE" 

16:  stp 

17:  beep; ent  "#  OF  FRAMES  TO  AVERAGE" , H 

18:  ent  "WINDOW  WIDTH  , FIRST  PIXEL" , r9 , "LAST  PIAEL",rlG 
19:  H+l-H 
20:  0-r3 

21:  for  0=1  to  H  by  4 
22:  for  1=2  to  9 
23 


35 

36 

37 
36 

39 

40 

41 

42 

43 

44 


24 

if 

1=2 

; "  930 

020" 

♦DS 

25 

if 

1=3 

"  020 

“  ♦O  $ 

26 

if 

1=4 

"93  0 

120" 

♦D$ 

27 

if 

1=5 

"  020 

"  ♦D  $ 

26 

if 

1=6 

"930 

220" 

♦D  $ 

29 

if 

1  =  7 

"  020 

"♦D$ 

30 

if 

1=8 

"930 

320" 

♦  D$ 

31 

if 

1=9 

"020 

"♦D$ 

32 

if 

len 

(D$) <4 ; gto 

38 

33 

char  (1] 

&"  2" 

&char 

(2) 

34 

cli 

7 

wait  100 
wrt  704,  BS 
cli  7 

cnar  (46) ♦C$ (1, 1] 

cnar( 1) &" 2  "&chac<2) &C$&cnar (15) &D$.  ( 1,3] 

cli  7 

wait  100 

wrt  704, AS 

cli  7 

if  I#3  or  0>l;gto  47 


■A$ 
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dsp  "BACKGROUND  RECORDED" 

wait  5U0;dsp  "CONTINUE  .vriEN  wX|  F RAHE  READY";  stp 
if  D$#  OlO"  or  1=3  or  1=9  ana  0<5;gto  49 
dsp  "CONT.  rtriEN  NXT  FRAME  SELEC TED " ; stp 

next  I 

if  ri-0+l<4;gto  53 
for  X=i  to  4 
gto  5  4 

for  X=1  to  H-G+l 

"  1  "&"  0"&"A"&"P"&cnar.  <96)  +  P$ 

if  X#l;gto  57 

1-»J  ;  0-*¥ ;  gto  62 

if  X»2;gto  59 

sfg  3;256-*J;  O¥;gto  62 

if  X#3;gto  61 

1*J ; 255-*¥ ; gto  62 

25o-*-J ; 255+Y 

for  L=1  to  256 

dsp  L 

0  +  D ;  "  "♦Q5;  Y  +  1-»Y 
Y*512+J-S 
for  1=2  to  6 
int(S/lo)  +T 

r-u 

3-T*lo+I*16*T 

U*S 

0$&cnar  (.X)  ♦0-$ 
next  I 

char  (1)  &"  =  "&cnar  (34)  &cnar (2)  &05&char  (3)  &P  5&char.  (15)  ♦  ES 
tint  1,  z,  c 
cli  7 

wrt  7U4.1,B5 

buf  "DATA" 

tfr  704, "DATA" ,256 

rds ( "DATA") *E; if  E=-l;jmpO 

for  I=r9  to  rlO 

num(  15  (1 , 1  ] )  *Z 

D+Z+D 

next  I 

D/  (rl0-r9) *R[L] 
next  L 

if  fig 3=0 ; gto  9 1 

"TAKE  DIFFERENCE  OF  IMAGE  AND  BACKGROUND" : 
rread  3,l;sread  3, 1(*J;for  1=1  to  256 
abs( R( I ] -I [ I ] ) ♦Rl I ] ; nex t  I 
gto  93 

rread  3,l;sprt  J,R[*J 
gto  113 
0*G ;G-F 

for  1=30  to  226 
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95:  R ( I  j  +G*G ; next  I 

96:  for  1*30  to  226  ;R(  I  ]  +  F*  F 

97:  if  F>*G/2  ;  gto  99 

98:  next  I 

99:  I* A 

100:  if  0*1  and  X*2;l*rl 

101:  rl-A*r2 

102:  for  J*1  to  25o+r2 

103:  R(  J-r2]  *R(  J J 

104:  if  J<256+r2; next  J 

105:  for  K*25o+r2  to  25b 

106:  0*R(KJ 

107:  next  K 

108:  sprt  7,r2 

109:  rread  4,l;sread  4,I(*] 

110:  for  1*1  to  256 

111:  I { I ] +R f I ] *R [I]  ; nex t  I 

112:  rread  4,1;  sprt  4,R[*J 

113:  next  X;next  0 

114:  "TAKE  AVERAGE  OF  ALL  FRAMES "  : 

115:  for  1*1  to  256; R( I ]/ (H-l) *R( I ] ;R( I ]*R(513-I ] ;nex t  I 
116:  rread  4,1; sprt  4,R(*J 
117:  if  0[15] <1; gsb  "OUTPUT" 

118:  beep 

119:  ent  “RANGE" , R  ,"SQCN" ,Q , "EXT  CCEF" ,£, "wAVLNGTH" , w, "OBJ  LE 
NS"  ,0 

120:  ent  "OBS/OBJ",  B, "SCALE" ,rll; sprt  6 ,R, Q, E , W  ,0  ,B, r 11 
121:  "SYSTEM  DATA-W*wAVLNTH  IN  METERS  O-DIA  OF  OBJ  LENS  IN  ME 
XERS" : 

122:  "B*DIA  OF  OBSUR/DIA  OF  OBJ  LENS  rll=SCALE  OF  DATA  IN  MIC 
RORADS/PT"  : 

123:  "R-DISTANCE  TO  rGT  IN  METERS"  : 

124:  "J-SQUARE  OF  ATMOS  TURB  INDEX  IN  METERS'* -2/3 "  : 

125:  "E-EXTINCTION  COEFF  (1/METERS)": 

126:  if  O [ 1 ] # 5 ; g to  140 
127:  beep 

128:  ent  "AMPLITUDE", A, "SIGMA”, C 

129:  for  1*1  to  256 

130:  rll*  ( 1-1) ♦F 

131:  F'*2/  (2*C'*2)  -*G 

132:  if  G>13 ;gto  134 

133:  A*exp  (-G)  ♦I  II  J ;  I*K;  gto  135 

134:  0*1(1] 

135:  next  I 

136:  rread  5,l;sprt  5,I(*];rread  5,l;sread  5 ,  R  [  *  J 
137:  1*0(15]  ;if  O.(14]#0;gsb  "OUTPUT" 

138:  gto  149 

139:  "SOURCE  PATTERN  DATA  INPUT": 

140:  rread  2,l;sread  2,R(*] 

141:  for  1*1  to  256 
142:  I*K 


143:  if  R( I ) <-500; gto  146 

144:  next  I 

145:  gto  149 

146:  for  J=I  to  256 

147:  0+R(J] 

148:  next  J 
149:  0(l]*i 

15U:  "BRANCH  ACCORDING  TO  TYPE  OF  SOURCE  PATTERN  DATA": 

151:  gto  152 ; if  I>l;gto  155;if  I>2;gto  161;if  I>3;gto  206;if 
I>4;gto  152 

152:  "CONVERTS  SOURCE  POINT  SPREAD  FUNCTION  TO  DINE  SPREAD  FC 
N"  : 

153:  gsb  "LSF" 

154:  2+0 { 15 J ; if  O.l3]>0;gsb  "OUTPUT" 

155:  "CALCULATION  OF  FOURIER  TRANSFORM  OF  LINE  SPREAD  FCN  OF 
SOURCE" : 

156:  for  1=1  to  256 
157:  R( I ] +R(513-I]  ;next  I 
153:  gsb  "FXFORM" 

159:  rread  5,l;sprt  5,R{*J 

160:  3+0  [15];  if  O[4J>0;gsb  "OUTPUT" 

161:  if  O[2]=2;gto  173 

162:  "OPTICS  DIFFRACTION  INPUT  IF  NOT  CALCULATED" : 

163:  rread  8,l;sread  8,R[*] 

164:  for  1=1  to  256 
165:  I+K 

166:  if  R[I ] <-500 ; gto  169 

167:  next  I 

168:  gto  182 

169:  for  J=I  to  256 

170:  0+R(J];next  J 

171:  gto  182 

172:  "CALCULATION  OF  DIFFRACTION  LIMIT  POINT  SPREAD  FCN": 

173:  B*2+0 ; 1-D+H; 3.14159e-6* rll*0/w+Z 

174:  for  1=1  to  256 

175:  Z*  (1-1)  *5f ; if  Y>30;gto  178 

176:  (  'AIR¥'(*)-D*'AIRK'U*B)  )  “2/d‘*2*Rl  IJ 

177:  I+K; gto  179 

173:  0+RUJ 

179:  next  I 

160:  rread  8,l;sprt  8,R[*J 

161:  4+0 [15] ; if  O.[6J>0;gsb  "OUTPUT" 

182:  if  0{2]=l;gto  138 

183:  "CONVERTS  OPTICS  POINT  SPREAD  FCN  TO  LINE  SPREAD  FCN": 
184:  rread  8,l;sread  3,R[*] 

135:  gsb  "LSF" 

136:  rread  8,l;sprt  3,R[*j 

1*7:  5+0 [  15]  ;if  0[7]>0;gsb  "OUTPUT" 

188:  "TAKES  FOURIER  TRANSFORM  TO  GET  OPTICS  MODULATION  TRANS F 
£R  FCN": 
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Ib9 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

21o 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 
227 
223 
229 


rread  3,l;sread  8,R(*] 

for  1=1  to  25o 

R[I J -r [ 513-1 J  ;nex t  I 

gsb  "FXFGRM" 

for  1=1  to  256 

if  R[l]=0;gto  196 

R  [  I  ]  /R  ( 1  ]  +R(  I  ]  ;next  I 

rread  8,l;sprt  o,R(*] 

6+0  [  15]  ;  if  0  [  8  j  >0;  gsb  "OUTPUT" 

"CALCULATION  CP  THE  PRODUCT  OF  TWO  FOURIER  TRANSFORMS" : 

rread  8,l;sread  8,R[*];rread  5,l;sreaa  5 , I [* } 

for  1=1  to  256 

R{ I ] *1 [ I ] +R[ I ] ; nex  t  I 

rread  10,1; sprt  10,R{*J 

7*0 [15]; if  O C 5 J >0; gsb  "OUTPUT" 

exp( -E*R) +X 

"ATMOS  MODULATION  TRANSFER  FCN  CALC-SHCRT  TERM" : 
0*R*21.6+A; W~ (-. 333333 3) *B ; G+V; 0+L; X*R[1] 
for  1=2  to  256 

gto  214;  if  L<  =0  ;  gto  215;  if  L=0;gto  20  9 

976.5625* (1-1) /rll*T;A*T~l. 666667* (S- (T/O) * .3333333) +G 

sfg  14 

if  013 ;  1*L 

if  G-V<0;-1+L 

G+V;X*exp(-G) +R(I] ;I*K;c£ g  14;gto  218 
0+R[I];gto  218 

R[K  J  *coa< (I-K) *3 . 14159/K) *2+Ri I j 
gto  217; if  I-3*K/2<0;gto  213 
I +K  ;  1 + L 
next  I 

rread  9,1; sprt  9,R[*] 

3+0  [  15  ] ;  if  O  1 9J  >0;  gsb  "OUTPUT" 

"TRANSFER  FCN  OF  SOURC£*OPTICS*ATMOSPHERE" : 
rread  10,l;sread  10, R]* J ;rread  9,l;sread  9,1 [ *  J 
for  1=1  to  256 
I[I]*R[I]-im 
I. [I]  *R[I] 

R[ I) +R{5I3-I] 
next  I 

9+0[  15]  ;  if  O  [  10]  >0;  gsb  "OUTPUT" 

"INVERSE  FOURIER  TRANSFORM  GIVES  TARGET  LINE  SPREAD  FCN" 


230;  sfg  l;gsb  "FXFORM"; I+K;cf g  1 
231:  10*O[15];if  Q[ll]>0;gsb  "OUTPUT" 

232;  "CONVERTS  LINE  SPREAD  FCN  TO  POINT  SPREAD  FCN  BY  ABEL  TR 
ANSFORM" : 

233:  gsb  "ABEL" 

234:  11+0[15]; if  Q(12]>G;gsb  "OUTPUT" 

235:  "CALCULATES  THE  FRACTION  OF  THE  POWER  INSIDE  A  CIRCLE  Or 
RADIUS  R" : 
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236:  gsb  "S’* 

237:  12*0{15];  if  Q(13]>0;qsb  "OUTPUT" 

238:  "READ  IN  NEW  RANGE  AND  Q  If  RANGE  NEG.  READ  IN  ALL  NEW  D 
AT  A"  : 

239: ""STARTING  WITH  I/O  CONTROL  CARD, IF  RANGE=G  PROGRAM  STOPS 
»  . 


240 
E  NL 

241 

242 

243 
24  4 
245 
24  6 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 
27  3 
27  4 

275 

276 

277 

278 

279 

280 
281 
282 
283 


"If  RANGE  IS  POS.  THE  PROGRAM  CARRIES  OUT  CALC'S  WITH  -TH 
»"  : 

"RANGE  AND  Q“ : 
if  O{l]=0;gto  245 

ent  "RANGE", R, "SwCN" ,0, “EXT  COEF" ,E 
gto  206; if  R<*0;gto  119;if  R*0;gto  245 
end 

"ABEL"  :R[1]*N;1.4R[1]-1.8R[ZJ  +  .4R(3]*RUJ  ;for  1=2  to  K 
R  [I  ] *M 

. 4N+. 2M- .6R [ 1+1  j  *R{ I ] ;M*N ;next  I 
for  1=1  to  K 

R]  I]  /  (2*  V(  (I  +  .1)“2-II))*R[IJ 
for  J=I  +  1  to  K 

R[  I]  +R[J]/V(  (J  +  .1)“2-II)*R[IJ 
next  J  ;R[  I] /ir*R[I  ]  ;dsp  I;next  I 
ret 

"LSF":for  1=1  to  K 
l*J;dsp  I ;R[ I] *Z 
|(I*I+J*J)  *Y 

2*  ( (1-f  rc(Y)  )*R[ int (Y)  3  +  frc  (Y)  *R[int  ( Y)  +1J )  +Z*Z 

J+l*J;if  Y<K ; jmp  -2 

Z*R[Ij;next  I; ret 

"B"  : .  2  5ffR  [  1]  *R(1] 

for  1=2  to  K 

2irIR[  I  ]  +R [  L-l]  *R[  I  ] 

next  I 

if  R[K]=0;ret 

for  1=1  to  K;R[I.]/R[KJ*R[I]  ;next  I;ret 

"AIRY":if  pl<0; beep;dsp  "ERROR  pl<Q";stp 

if  pl=0 ; l*r4; ret  r4 

0*r5;if  pl>15;jmp  2 

20+10pl-pl“ ( 2/3) *r6 ; jmp  2 

90+pl/2*r6 

if  pl<5;6+pl*rl2; jmp  2 

1. 4*pl+6Q/pl*rl2 

max ( int (r 12) ,int(3+pl/4) )*rl2 

for  M»rl2  to  r6  by  3 ; le-28*rd ;0*rl3*rl4 

sfg  10;if  M/2  =  int(M/2)  ;cf g  10 

for  J»1  to  M-2;2*(M-J)  *r8/pl-r  13*rl5  ;  r  8*rl3 

rl5*r3;if  M-J-2=0 ; rl5*r 4 

cmf  10 ; rl4+2*r8*f lgl0*r 14 ; next  J 

2*r8/pl-rl3*rl5 

rl4+rl5*rl4 ; r4/rl4*r4 

if  abs(  r4-r5) -abs (r 4*le-6)  <*0  ; 2*r4/pl*r4;  ret  r4 
r4*r5;next  M 
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beep;dsp  "Accuracy  not  obtained" ;wait  i000;ret  r4 
"FXFORM" :rad;9»N;ina  I 
n/2  *  (N-l )  »T 

for  M*1  to  N; 2" (N-M) *rl6 

for  J=0  to  2*{M-l)-l;cll  'a I ' ( J, P ,N-1) 

cos(PT)-C;sin(PT)  (l-2flgl)*P 

for  I*2r  16J  +  1  to  2rl6J  +  rl6 

R  [  I]  -*>r0 ;  R[  I+r  16 ]  ♦rl 

I[I]+r2;I{ I+rl6]*r3 

rO+rl*C+r3*P«-R[  I] ; r2+ r3 *C-rl*P-I [I] 

rO-rl*C-r3*P->R[I+rl6] 

r2-r3 *C+r 1*P+ I [ 1+ rl6] 

next  I;next  J;next  M 

for  1*0  to  2"N-1 ; ell  'BI  '  ( I ,  J,  N)  ;  if  I-J>0;jmp  8 

if  I*J;jmp  5 

R[I+l]/V(2“tf)-P 

1. 1 1  +  1]  /  V  (  2*N )  ♦Z 

R[J  +  1]*R[  1+1]  ;  I  [  J+l  ]  ♦L(I  +  1J 

P-R[J+1] ; Z-»I { J+1J 

R[ 1  +  1]/ V ( 2~N) -R[ 1+1] 

I[I+1]//(2*N)*I[I  +  13 
next  I;deg;ret 
"BI“  j0*p2;  pl-»p4 
for  Z*1  to  p3 
p4/2*p4;  2p2*p2 
if  frc(p4)#0;p2+l*p2 
int(p4) -*p4; next  Z;ret 
"SOURCE" :sf g  2 

beepjdsp  "DISPLAY  SOURCE  IMAGE  ON  VID  DISC"; wait  30( 

ell  'SCALE' 

rread  2,l;sprt  2,R[*] 

cfg  2; ret 

"SCALE" :if  f lg2; gto  318 

beep;dsp  "DISPLAY  SCALE  IMAGE  ON  VID  DISC"; wait  300( 
dsp  "THEN  PRESS  CONTINUE Stp 

ent  "WINDOW  WIDTH,  FIRST  PI XEL" , r9 , "LAST  PIXEL", rlO 
for  1*2  to  3 

char (47}&char ( 55) &char  (6  4) &char (8U)&char ( 96)  +  C$ 
if  1*2; "930  020 "-DS 
if  I*3;"020"+D$ 
if  len(D$)  <4  ;  jrap  5 

char  (1)&"  2"&char  (2)  &C$&char  (15)&D$(5,7]*B$ 

cli  7 

wait  100 

wrt  704, B$ 

char.  (46)*CS  (1, 1] 

char  (L)&"  2  "&char.  (2)  &C$&char.  (15)  &D$  (1  ,.3)  +A$ 
cli  7 
wait  100 
wrt  704, A$ 


DISC" ; wait  3000 


VID  DISC"; wait 


334 

335 

336 
3  37 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 
.353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 
.368 

369 

370 

371 

372 

373 

374 

375 

376 

377 
373 

379 

380 

381 

382 
363 


next  I 

"  1  0"  &  "  A"  &  "P"  &  char-  (9  6)  »P$ 

1*J ;0*Y 

for  L=1  to  256 
dsp  L 

0*D;  ""♦OS;*  +  l-** 

**512+J*S 
for  1=2  to  6 

in t (S/16) *T;T*U;S-T*16+I*lo*T 
U-*S;G$&cnar  (.T)  ♦0$ 
next  I 

char  (L)&"»"&char  (34)  &cnar  (2)  &0$&char  (  3)  &P$&char  (15)*B$ 
cli  7 

wrt  704.1, B$ 

buf  "DATA” 

tfr  704  ,  "DATA", 256 

rds  ( "DATA")  ♦£ ; if  E»-l;jmp  0 

for  I=r9  to  rlO 

num( I  $ (I ,1 ] )*Z;D+Z»D 

next  1 

D/(rlO-r9)*R[L] 
next  L 

if  flg2;gto  358 

if  0[16]  =2;  gsb  "OUTPUT" 

ret 

"OUTPUT"  : 
beep 

dsp  "STOP  IF  PLOTTER  NOT  READ*  THEN  CONTINUE” ;wait  5000 

sfg  14;''MIN'(256)/'riAX'{256)-rl7;cfg  14 

if  rl7<-.l;scl  75 , 1. 05,-30 , 264 ; pit  l,0,l;jmp  2 

sc 1  -.15, 1.05, -30, 264; pit  1,0,1 

for  1=10  to  0  by  -1 

pit  1/10,0, 2;plt  1/10,256/150,2 

pit  1/10 ,0,2; next  I 

for  1=0  to  16 

pit  0 , 161 , 2; pit  1/150,161,2 
pit  0,161  ,2 
next  I 

pen;csiz  .9, 1,1. 4 

fxd  0;for  1=16  to  0  by  -l;plt  -.1,161-1,1 
lbl  16(16-1) ;next  I 

fxd  l;for  1=0  to  10;plt  1/10- . 06 ,-5 , 1 
lbl  I/10;next  I 

if  O [15 ] =10  or  O [ 1 5 ] * 0 ; " TARG ET  LINE  NUMBER"+*$ 
if  O  [15]  =1  or  Q[  15]  =11;  "POSITION  VALUE"+*$ 
if  0(15] =10 ; "PREDICTED  NORMALIZED  LINE  VALUE"*X$ 
if  0(15] =0; "MEASURED  NORMALIZED  LINE  VAL UE"*X$ 
pit  -.12,51, l;csiz  1. 2 ,1 , 1. 4 , 90 ; lbl  *$ 
pit  0,-13,l;csiz  1. 2 , 1, 1. 4 ; lbl  X$ 
csiz  .5,1, .7,0 
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334:  'MAX'(25b) *r7 

385:  if  r7=0;l*r7 

386:  for  J=1  to  256 

387:  pit  ft( J]/r7 , 257-J;nex t  J 

388:  penjret 

389:  "MIN"  :R[l]--p2 

390:  for  1=2  to  pi ; if  R[I]>p2;jmp  2 

391:  R ( I ] ♦p2 

392:  next  I; ret  p2 

393:  "MAX"  :RUJ*p2 

394:  for  1*2  to  pi; if  R[I]<p2;jmp  2 

395:  R[I]  *pr. 

396:  next  I;ret  o2 
*21194 
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